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The Anderson-Holstein model is studied in the framework of the semiclassical approximation. 
Analytic results for Kondo temperature renormalized by weak electron-phonon interaction and for 
phonon Green function are obtained, and they are interpreted from the viewpoint of dynamical mean 
field theory. Especially the isotope effect of the effective electron mass is discussed in the presence of 
strong electron correlation. The results are also compared with those by numerical renormalization 
group and other related works, and they are consistent with each other in their common domain of 
validity. 
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I. INTRODUCTION 

Both electron correlation and lattice dynamics are im- 
portant in understanding the physical properties of many 
condensed matter systemsi But the interplay between 
them has not been studied in detail owing to the lack 
of reliable calculation methods. Conceptually this prob- 
lem is difficult to handle due to the absence of suitable 
expansion parameters. 

Impurity problems have been relatively well studied 
compared to problems on lattice, mainly thanks to di- 
verse analytic and numerical nonperturbative techniques. 
With the recent advent of the dynamical mean field the- 
ory (DMFT)2 various lattice problems can be mapped 
(or approximated) to impurity problems. The mapping 
can be justified rigorously in the limit of infinite spa- 
tial dimensions. Once the mapping to a certain impu- 
rity problem is accomplished the powerful techniques of 
impurity problems can be employed to understand the 
lattice problems. 

In the perspective of DMFT the problem of the inter- 
play between electron correlation and lattice dynamics 
can be mapped to an impurity problem with both elec- 
tron and phonon degrees of freedom. One of the simplest 
model with both degrees of freedom is the Anderson- 
Holstein (AH) models AH model is a single-impurity 
Anderson model with a linear coupling to a local phonon 
mode as in the Holstein models 

Our goal is to understand the physical properties of 
the AH model in the light of the interplay of electron 
correlation and lattice dynamics. Recently the numerical 
renormalization group (NRG) method has been success- 
fully applied to AH mode&4, and almost exact results 
on the electron and phonon spectral functions have been 
obtained. However, it is still desirable to develop an an- 
alytic scheme in spite of its approximate nature since it 
helps us to understand the underlying physics in more 
clear and intuitive way. Of course, a systematic analytic 
approach to a model like AH model does not exist in gen- 



eral. But by restricting our focus to a certain parameter 
regime we can develop an approximation scheme which 
captures the essential features of the interplay between 
electron correlation and phonons. We will be mostly in- 
terested in the regime where the electron correlation ef- 
fects are much stronger [large local Coulomb repulsion 
U] than the electron-phonon coupling. If the electron- 
phonon coupling is omitted (Anderson model) it is well 
known that the ground state is the Kondo singlet state. 
Our main focus is on understanding how the spin fluctu- 
ation processes responsible for the Kondo ground state 
are influenced by the weak electron-phonon interaction. 
We have developed a semiclassical approximation scheme 
in the context of path integral approach where the spin 
fluctuations are represented by instantons ( or kinks). 6 In 
this scheme phonons can be naturally incorporated and 
the features of interplay can be clearly exhibited. 

Isotope effect is a suitable physical property for ad- 
dressing the issue of the interplay of electron correlation 
and phonons. In conventional metals, where the Midal- 
Eliashberg (ME) theorjii of electron-phonon applies, the 
isotope effect appears to be very small. This is because 
ME theory keeps only the diagrams of the zeroth order 
of Q/t, where fi is phonon frequency and t is the band- 
width of conduction electron^ The omission of diagrams 
of higher order in 0,/t is justified by the Migdal theo- 
rem on the electron-phonon vertex function, which is the 
key element of ME theory. Presumably the ME theory 
cannot be employed for the systems with strong electron 
correlations since vertex corrections are not expected to 
be small for such systems in general. We can compute 
the isotope effect for the strongly correlated systems by 
studying how the Kondo temperature of impurity model 
is modified by (local) electron-phonon interaction [see Eq. 
©]■ 

There is a caveat in directly translating the results ob- 
tained for impurity model to those of the corresponding 
lattice model. The lattice model with Lorentzian density 
of state (DoS) is formally identical with the impurity 
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model since there is no need for solving self-consistency 
relation. Thus for the lattice model with the Lorentzian 
DoS the results obtained for the impurity model can be 
directly applied. However, the Lorentzian DoS is not 
physical in many aspects, but it is expected to be valid 
for the systems with finite DoS at Fermi surfaced The 
constraint of finite DoS at Fermi surface restricts the 
applicability of our results to the metallic regime. For 
the lattice Hubbard model this implies that the system 
should be away from half-filling or [/ c2 > U > 7rA at 
half-filling [see Ref. .2] for details]. In addition, our re- 
sults are not relevant for the insulating regime such as 
small polaron regime^ 

We point out that Deppeler and Millis addressed the 
similar issue using the adiabatic expansion of DMFT 
without assuming any specific ground state&ifi They 
identified a new class of electron-phonon diagrams missed 
in previous studies but they treated electron correlations 
in the scheme of Hatree-Fock approximation. Our ap- 
proach treats the strong electron correlation more accu- 
rately, thus it provides the more accuruate electron cor- 
relation dependences of physical properties. 

In a single band model without extra orbital de- 
grees of freedom, the Holstein phonon interaction acts 
in the charge channel, while the electron correlation ef- 
fect, represented by local Coulomb repulsion, acts via 
the spin channel at low energy. Due to this apparent 
"spin-charge separation" , the interplay of electron cor- 
relation and electron-phonon interaction first appears in 
the second order expansion with respect to the phonon 
field in our treatment. For a multiband model where 
the orbital degrees of freedom is active electron-phonon 
dynamics becomes much richer as exemplified by the 
phase diagrams of colossal magnetoresistance (CMR) 
materialsiii We have compared our results with those of 
the spin-polarized Jahn- Teller model studied in the con- 
text of CMR manganites^ The only formal difference 
between two models lies in the structure of the coupling 
of phonon modes, but this difference brings about drastic 
differences in physical properties such as isotope effects. 

Our main results are the following: (1) the electron- 
phonon interaction increases the Kondo energy scale in 
the regime of strong electron correlation. The explicit 
results are given in Eqs. (|31)1I3§|) . (2) the isotope effect of 
effective electron mass Eq. Q44[l. (3) the phonon spectral 
function computed in the fugacity expansion scheme, Eq. 

This paper is organized in the following way: In Sec. II 
we introduce AH model and present it in a form suitable 
for semiclassical approximation. In Sec. Ill we develop 
an instanton expansion of partition function and discuss 
its renormalization group (RG) flow. In Sec. IV the 
phonon Green function is calculated in the fugacity ex- 
pansion scheme. We conclude this paper in Sec. V with 
discussion and summary. Some details of calculations can 
be found in appendices. 



II. FORMULATION 



The AH model is defined by the following 
Hamiltonian^ 



'Hah — Ti-ei + Hph + 7~Lei- P h 



(1) 



H ph = \(MtfQ> + ^). 



H e i- ph = gQC£flf < r-l). 



(3) 



(4) 



a =1,1 is the spin index, and V k is the hybridization 
matrix element. f a is the impurity electron operator, 
and Cfc CT is the conduction electron operator. Q is local 
phonon coordinate. P is the conjugate momentum of Q 
satisfying [Q, P] = ih. M is ion mass and SI is the oscilla- 
tor frequency of dispersionless (Einstein) phonons. Ai at 
is the number of lattice sites for the conduction electrons. 
We will consider symmetric AH model specified by the 
relation e f + U/2 = 0. 

The conduction electrons c& can be integrated out ex- 
actly, yielding local impurity action. 



r/3 

S imp = -TV fl G 1 (ie) f a + dr Uf\frf\ 

Jo 



fl 



dr 



Q 2 +gQ(n-l) 



(5) 



where n a = n = n T + ■ 



and 



G \ie) = ie-e/-S(«e), 



E(ie) 



— r 

AT,.. ^ 



Mat ^ ie + fi- e k 



(6) 



For the hybridization matrix elements weakly dependent 
on momentum near Fermi surface, G~~ l (ie) can be ap- 
proximated by 



G^{ie) = ie- e f + iAsgn(e), 
where A is the hybridization enery scale 
A = ^(H4| 2 ^-e fe )) avg . 



(7) 



(8) 



The average is done on the Fermi surface of conduction 
electron. We are interested in the strongly correlated 
regime U ^> A. The charge fluctuations at energy scale 
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of order U 3> A can be neglected in considering low en- 
ergy physics. To separate out the charge fluctuations at 
high energy it is convenient to introduce the Hubbard- 
Stratonovich (HS) transformation for the Hubbard U in- 
teraction. Two HS fields <fi c and <p s are introduced to 
decouple on-site repulsion U into charge and spin chan- 
nel, respectively. 



e -un m _ e - T 



D\ 



c ,<ps\ exp 



U 



V 



+ icf> c n + 4> S S Z ]. (9) 



S z = nj — ni is the z-component of impurity spin. Then, 
we choose the mean value of charge HS fields which 
can be absorbed into G^" 1 (ie) [cancelling — e/] and ne- 
glect its high energy fluctuations. 6 Then electron-electron 
interaction [valid only at low energy] acts only in the 
spin channel, while phonon couples only in the charge 
channel. Consequently at the leading approximation 
electron-electron interaction and electron-phonon inter- 
action seem to decouple, however the interplay can be 
seen from the second order in expansion with respect to 
phonon fields [see below]. 

The resulting local action in imaginary time is 

Simp = - J drdr'J^fUr^GoHr' -r)f a (r) 



dr 



r 1 

U 



4> 2 S - s (n T - ni) + 5<5(n T + n l - 1) 



M 



-(d T Q) 



, Nin 2 , 

Q 



2 v . -w 2 
The bare Green function Gq(t' — r) is given by 

g — «u)„(t'— t) 



G (t'-t)=t£ 



iuj n + zAsgnu;,, 



(10) 



(11) 



The impurity electron / can be integrated out exactly 6 
for the action of the type of Eq. (fTU|) . 



D[f, <f>.,Q]e 
where S = Sp + Sk is given by 

^ r^2 M 



D[<f> s ,Q}e- s ^^\ (12) 



u 



-n z Q 



£ ^ ( - t* tan" 1 & + \ ln[l + £])] , (13) 



convenient to define an energy scale El associated with 
lattice relaxation (or equivalently polaron energy). 



E L = 



2Mfl 2 ' 



(15) 



Following Ref. we define dimcnsionless variables: 

Y 4>s v gQ u e l n 

X = — . Y — — —, u——, A——— 7=—. 16 
A' A A A A 

To develop semiclassical approximation we first need to 
locate the classical minima of potential term Sp. When 
u 3> A and u > tt, Sp has two degenerate minima 
(X*,Y*) = (±770,0), where 770 ~ u/2. We also assume 
7 < 1, even though the adiabatic expansion is not sys- 
tematically employed in our study. Our main interest 
resides in the case with strong electron correlation and 
weak electron-phonon interaction. Accordingly, we re- 
strict our attention to the regime u 1 3> A, so that the 
phonon part can be treated perturbatively. 

There can exist other minima in the regime A 3> u. 
This regime is essentially similar to the original Holstcin 
model and has been studied in Ref. 01 m l ne semiclas- 
sical approximation scheme. Since we have treated the 
charge HS field <j> c , which couples to charge density like 
phonons, only in a mean-field approximation, our present 
approach does not apply to the regime A»u. Further- 
more, the polaron state which is likely to occur in the 
regime X^S> u cannot be described within DMFT scheme 
with Lorcntzian DoS. 

The quantum tunnelings between two minima 
(X*.Y*) — (±770,0) are nothing but the spin flip pro- 
cesses in the path integral framework The tunneling 
amplitude is determined by the kinetic term Sk- The 
key is to understand how the presence of weak electron- 
phonon interaction influences the spin flip processes 
which are responsible for strongly correlated states. To 
this end we integrate out the phonon field Y in the second 
order. There are no parts linear in Y. They vanish upon 
the spin summation. This integration over phonon field 
is justified since the phonon field has a trivial minimum 
at origin, so that we can employ harmonic approximation 
at low energy. Concretely, we have to expand the action 
Eqs. i|13l EH as a power series in Y up to the seoend 
order. The action S = Sp + Sk expanded in Y can be 
written as 



S ~ S x ± S' 



(17) 



f P , M ro v-^ 1 ft 3 drdr' 

LJ-to i+m (14) 

where Cer( r ) = [ — cr 0s( r ) + .9 ( 3( r )]/A. Sp and Sk is 
the potential and the kinetic term, respectively. It is 



Sx is just the action Eqs. (| 1 31 114[) with Q = 0, and it is 

the starting point of the Hamann's original anaylsis£ 



&x — — 



drdr 



tt Jo t — t 



7 X(r) 



111 1+X 2 (t') 



dX{r') 
dr' X 2 (t)~X 2 (t') 



dr 



\X 



1 



— + -{-X tan" 1 X + - ln[l + X 2 ]lj8) 
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Sy is the phonon part of Eq. I|10|) written in terms of Y . 



Sy = dr 
Jo 



4AA 7 2 I dr J 4A { ' 



(19) 



Sxy embodies the interplay of spin fluctuation and the 
electron-phonon interaction, and has very complicated 
structure. The explicit expressions of the terms of Sxy 
are listed in Appendix A. A typical term of Sxy looks 
as follows: 



S 



1 



XY,2 



drdr' \n\r -t'\Y{t)Y{t') 



/ drdr' In |t — r'| Y(t)Y(t') 



2£A 2 (t) 



+ £X 2 (t)) 2 (1 + £A 2 (t')) 



.(20) 



III. RENORMALIZATION GROUP ANALYSIS 

The partition function can be computed by summing 
over instanton trajectories [ or hopping paths] in semi- 
classical approximation^ The instanton trajectories are 
characterized by the hopping locations {ti}. For simplic- 
ity we choose the same trajectories as those of Hamann's, 
so that the derivative of X is strictly zero unless ti < r < 
ti + to. to is the range in imaginary time where the hop- 
ping occurs. We then expand the partition function Z 
in the number of hoppings N = 0, 2, 4, • • • . The num- 
ber of the hoppings is even due to the periodic boundary 
condition in imaginary time. 



Z — Z N=0 + Z N=2 + Z N=4 + 



(21) 



It is convenient to normalize the partition function with 
respect to Z^ =0 = Z a . Z a is the partition function for 
the classical minima X — ±770 • Since the action is an 
even functional in X one can choose either minimum in 
computing Z . We choose X = +7/o- 



Z = J D[Y] 



-S ph [Y]-S x [X=r, ] 



S x [X = ri ] = P 

l a 

ln(l + tjg) 



„2 2 

— + -(-770 tan 1 m 

U IT 



S ph [Y] = S Y + S XY [X = Vo ] 



i i 



AA7 2 



7r(l + ry 2 ) 2 ' 

)l|yM| 2 , (22) 



A 



41 + V 2 o)' 



Note two new terms coming from Sxy [A" = ?7o] i n 5ph [Y] . 
From the action S p h in Eq. (12211 the phonon propagator 



in the absence of instantons is easily obtained. 

M _1 



(Y(iu;)Y(-iu;)) N = = ) 



A' M 2 + £*M + cj 2 ' (23) 



where [ the condition 770 3> 1 is used] 

-0 2 - ^ 2 (i- A). 

4 7 2 

E ~ El t *C El, 



a = E*/A. 



(24) 



Note that the frequency of phonon decreases as the elec- 
tron phonon coupling g increases in the strong coupling 
limit of large U. This agrees with the result of Ref.Q. 
E* is a new energy scale involving both electron-phonon 
and electron-electron interaction. Note that E*/Cl — 
A 7 /^«l. 

Being expanded in the number of instantons the nor- 
malized partition function can be expressed as 



Z 
Z~o 



\2n 



dti n f t2n T ° dt2n-l 



T Jo 



l + ^T(zT 
n=l 

* t2 ~ T ° dti 

I — exp[+V(ti,-- ■ ,t 2n 







TO 



(25) 



where ^({tj}) is the energy of instanton gas which in- 
cludes the corrections induced by phonon coupling, z 
is the fugacity of instanton gas. Zq provides a denom- 
inator for phonon averaging, so that the explicit form 
of (Y (iu>)Y (—iu>)) n=o can be used in the evaluation of 
partition function. The integration over phonon field Y 
induces a few corrections to the original Coulomb (kink) 
gas action&i^of Kondo problem without phonons [K% 
term of Eq. 1|26|) ] . These corrections represent the inter- 
play of electron correlation and electron-phonon interac- 
tion. 



i>3 



ti tj 

TO 



+ ^ 2 ]Tln(^-^) d(t;. /,).(.+ ! 



i>j 



TO 



1 dY 

d(t t ,t 3 ) = -j (— (t = U) — {t = tj)). (27) 



dY 
~dr~ 



The intial (bare) values of coupling constants of Eq. i|26|) 

are 



K x ~ 2 
K 2 ~ - 



2 tan 1 770 
(^o/A) 2 



A 7 
To 



ztq ~ exp 



In 77 



A7 
2^ 



(28) 
(29) 
(30) 
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ci is a numerical constant of order unity and 

„2 

7 = w / A, A 



,9 Z 



2Mw§ A ' 



(31) 



The corrections in Eqs. (|28I3L)[1 [ Al ] reveal the expan- 
sion parameters as noted by Deppeler and Millis. 10 The 
influences of phonons are reflected in the K 2 term of the 
energy U({ti}) and in the change of initial values from 
those without phonons. The most drastic difference of 
K 2 term from K\ is the absence of the alternating factor 
(— , This absence of the alternating factor indirectly 
demonstrate the fact that the phonon interacts with elec- 
tron via charge channel not spin channel. 

The renormalization group equation (RGE) of the in- 
stanton gas Eq. (I25|l can be derived by following the 
procedure of Ref. |l4|. For details, see appendix II. 



(IX r, „ „ -1 

# = x[(l-^-) + 2x 2 K 2 g 2 ], 



dl 



dl 



(32) 
(33) 



where dl = dro/ro, x — zt^. g 2 is defined in Eq. l|B10fl . 
At the present level of approximation the influence of 
phonon coupling is manifest only in the renormalization 
of the fugacity x. For the sake of the comparison with 
the original Kondo model it is convenient to introduce a 
new variable y — 2 — K\ . The initial value of y is smaller 
than 1 as can be seen Eq. H28|) . Using the result Eq. 
(|B14j) . the RG Eqs. I|32I33|) can be re-expressed as 

dx x ' 

t: v 



dl 

dy 
dl 



x 2 ae 



4yx 2 ~ 8a; 2 



dK 2 
dl 



0. 



(34) 



a is defined in Eq. (|24|l . Since the phonons have been 
treated perturbatively from the outset we are not far from 
the original Kondo physics which does not have phonons. 
Thus, even if the SU(2) spin symmetry is not mani- 
fest in our approximate RGE, we can trace the isotropic 
ray y = Ax of Eq. I|34l) in approaching strong coupling 
regime. Due to the presence of explicit energy scale ojq 
in Eq. (|3lJ) . two energy scales Tjp , uj should be taken 
into account in the analysis of Eq. JSIJl. T^ 0) ~ De~ u / A 
is the Kondo temperature in the absence of phonons. D 
is an energy cutoff. 

First consider the case of luq <C TjP . Then the scaling 
should stop near Ik ~ In < In ^ = l Uo . This is 

because beyond Ik the RG flow begins to enter the strong 
coupling regime where our perturbative scaling equations 
do not apply. In this region of I < Ik, e~ u<)T ° ~ 1 holds. 
Integrating RG equation until x becomes order of unity 
we obtain 

T K ~ Dexp[-—^-— ] x exp[^ ln( 



2x(0) J 

X 1 °1°(iq) 



^(0)' 



(35) 



According to the result Eq. (|35|l the electron-phonon 
interaction increases the Kondo energy scale, and this 
is in agreement with the results of NRG studiesiii The 
amount of increase is given by 



ST, 



A 



A" 



(^2 + 8 Inryo)^ 

(0) 



(0) 



(36) 



Second consider the case of T K < wq. In this case, the 



RG flow evolves until I = L 



br- 



and then the sec- 



ond term of ^| in Eq. I|34|) becomes completely negligible 
beyond l U(j . In the region / > l Ua , x(l Uo ) plays a role of a 
new initial value. Integration RG equation starting from 
I = lui until x becomes the order of unity we obtain the 
Kondo temperature for the second case. 



T, 



K 



Woe -l/2x(^ 



(37) 



The direct evaluation of x(l UQ ) gives 



K 



T. 



(0) 



A" 



exp 



A a I 
—. j + — In — 

2% 2t] 70 



70 = uq/A. (38) 



The increase of Kondo temperature is given by 



ST, 



A 



A a , 1 
— 7 H In — 

2 Vo 2r )o 7o 



T, 



(0) 



A 



(39) 



Note that the results Eqs. I|35I38|) are obtained in the 
regime of strong electron correlation where ME theory 
does not apply. Both results Eqs. (|35I38|I show the 
electron-phonon interaction increases the Kondo temper- 
ature but by rather small amount. The increase of Kondo 
temperture could be heuristically understood in some- 
what unphysical antiadiabatic limit ^ A [see Eq. 
I62|l ] . This point will be elaborated further in sec V. In the 
antiadiabatic limit the retardation effects of phonons are 
negligible. Our results of Kondo temperature are valid in 
the strong electron correlation and fully incorporate the 
retardation effects of phonons. 

The results, Eqs. I|35I38|I . can be interpreted in the 
light of lattice system within the context of DMFT. The 
Kondo singlet state of impurity problem corresponds to 
the Fermi liquid in the original (large- d) lattice system 
(see the Sec. VII of Ref. 0). Here the momentum de- 
pendence of the electron self-energy is neglected, and the 
effect of correlation is encapsulated in the wave-function 
renormalization Z: 



Z{iu) 



1 



<9£(z) 



dz 



(40) 



which gives the residue of the quasiparticle. Z is related 
to the mass enhancements of quasiparticles as 



m* /mi,, 



(41) 



where m* is the effective mass, and rrib is the bare mass 
given by band theory. The Kondo temperature itself can 
be identified with the effective bandwidth B Q g of lattice 
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system. Then the mass enhancement can be estimated 
to be 



Z~ l = m*/m b ~ D/T K . 



(42) 



D is the order of half bandwidth. In the DMFT with 
Lorentzian Dos it can be identified with A. The above es- 
timate follows from the observation B e g ~ p F /m* , D ~ 
p 2 F /m. pp is Fermi momentum. Combining Eq. I)42|) 
with our results Eqs. i|35l38l) . the issue of isotope effect 
on the mass enhancement can be addressed. The ratio 
of mass enhancements in the presence and absence of 
phonons is given by 



m*/mb 
ji*>(°) /nib 



TP/T K ~ 1 - ST K /T } 



(0) 
K ' 



(43) 



where tyi* * is the effective mass in the absence of 
phonons. Thus, the isotope effect is negative in that 
the electron-phonon interaction which depends on ion 
masses decreases the effective mass of electrons. This 
isotope effect manifests itself in the specific heat coeffi- 
cient, the T 2 coefficient of low temperature resistivity, 
the Drude weight of optical conductivity, and in other 
physical quantities which depend on the effective band- 
with B b r in a crucial way. Our results show that the 
isotope effect is very small in the regime of strong elec- 
tron correlation. 



m* /mi, 
n*'(°> /nib 



1 - E L A/U 2 



(44) 



The proof of Eq. "4~4l is beyond the validity of ME theory, 
however. 



IV. PHONON CORRELATION FUNCTION 



The phonon correlation function is defined by 

D(t _ t>) = lD[X,y}e- S Y(t)Y( t ') 
1 ' jD[X,Y] e s ■ 

The correlation function of original phonon coordinate Q 
is related with Eq. (|45|l by 



(46) 



(Q(t)Q(t')) = (-fD(t-t'). 
9 



The denominator of the right hand side of Eq. (|45|l is 
just the partition function stuied in the previous section 
and will be written as Z = Zq(Z/Zq). 

We will compute the phonon correlation function in 
the framework of fugacity expansion. The phonon cor- 
relation function in the absence of instanton is given by 
Eq. l(2*37 . As nas been discussed in Sec. HHI the fugacity 
is relevant at low energy, thus our result for the phonon 
correlation is expected to be valid only in the high fre- 
quency/temperature regime (|w|,T) > Tk- 



The first nontrivial contribution comes from the trajec- 
tories with an instanton-antiinstanton pair. We rewrite 
Eq. (gSJ as 



D(t - = 



z ph I D[X,Y]e- s+s ^Y(t)Y(t') 



(47) 



where Z p h = J D[Y]e~ s ^ h[Y] [see Eq. (j22J]- Let us write 
the exponent of Eq. H47|) in the form 



-S + S x [vo] = -S p h - SS_ 



x 



SS 



XY, 



(48) 



where 



SS X = S x [X\-Sx[X = T)o], 
SSxy = S X y[X,Y]-Sxy[X = t)o,Y]. (49) 

By construction SSx and SSxy vanish in the absence of 
instantons. Now the structure of fugacity expansion is 
manifest in the form of Eq. (|47|l . 



D(t - t') = D (t - t') + D 2 {t - t') + 



(50) 



where D n is proportional to z n . Since we are only in- 
terested in the connected diagrams we have to select the 
connected components among contributions to Di(t — t'). 

The structure of SSxy is very complicated, and a typ- 
ical term ( which, in fact, will give the most dominant 
contribution ) looks like 



SSxy,2 = — - 



1 



1 



drdr' \tl\t -t'\Y{t)Y{t') 



1 



i 



1 + £X2( T )1 + £XV) (i+^o 2 ) 2 

1 



+ 4t / drdr' In It - t'\Y(t)Y(t') 

7T 2 



2£* 2 (r) 



d£ 







(l + ^ 2 (r)) 2 (l + ^ 2 (r')) (l + e% 2 ) 3 



(51) 



Next we need to evaluate the action SSx and SSxy in 
the subspace of two instantons which are parameterized 
by their locations t%,t2- Then the connected component 
of Diij — t') can be written 



D 2 {t-t') 



P h 



D\Y]e" 



rtz—Ta 

dti I dti 

n JO 



x e - Klln — e- SSxY ^^Y(t)Y(t'). (52) 

We will treat 5SxY(ti,t2) perturbatively as in the eval- 
uation of the partition function. In the expansion 
e -ss X Y — 1 _ SSxy + ■ • • , the first term cancels the 
disconnected piece coming from the denominator Z/Zq. 
Among various contributions to D 2 (t — t') the term con- 
taining two time derivatives of Y, 5Sxy,2, gives the most 
dominant one. The calculation of D^it — t') can be cast 
into the following form 

D 2 (t - t') = D (t -t') + J dhdt 2 
xD {t-t 1 )V{t 1 -t 2 )D {t 2 -t'). (53) 
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The function £(ti — £2) plays a role of self-energy. The 
explicit form of self-energy is 



E(t 1 -t 2 ) - ^-r 



^4 / Iti-faK jf, 

V Til ' 



(54) 



In frequency space it becomes 

,2 



EM-z^M^ln^ 
% M 



(55) 



Thus the phonon correlation function which is accurate 
up to the second order of z 2 is 



D(iw) 



2A 7 2 A 



ui 2 + E* \w\ + w 2 + z 2 i 1 E 



Apparently the self-energy term behaves very singular 
at low energy, but we have to keep in mind that the 
Eq. (|56[l was derived in the high temperature/frequency 
regime (|w|,T) > T K . 

Let us compare our results with those by NRG. The 
NRG calculation of phonon Green function by Hewson 
and Meyer has been improved in Ref.^jJ. The imaginary 
part of phonon Green function in the regime of large U 
has the following features: (1) the decrease of phonon fre- 
quency with increasing phonon coupling (2) broadended 
phonon peak (3) absence of peak near zero frequency. 
These feature are consistent with our results. First the 
renormalized phonon frequency 



U! 



(57) 



clearly shows the decrease of phonon frequncey with in- 
creasing A. 

This result can be compared with that derived in the 
weak coupling RPA type perturbative calculation^ 



UJ 



n i - 



El 

7rA 



(58) 



Second the renormalization of the phonon propagator 
can be reliably discussed if the interested energy scale is 
larger than Kondo temperature. If uj > Tx our result for 
the phonon Green's function can be applied in the neigh- 
borhood of wo. The broadening is basically determined 
by the energy scale E* which increases with increasing 
electron-phonon coupling. This feature is also consistent 
with the NRG result. In the perturbative RPA approach 
the broadening width is predicted to be linear in u> while 
our results predict sublincar dependence of width uj 1 ~ y 
which is stronger at low energy. Notice that the anoma- 
lous self-energy term represents strongly temperature de- 
pendent broadening since it originates from the electronic 
processes happening near Fermi energy. This feature can 
be checked in principle by NRG calculation at finite tem- 
perature. 



V. DISCUSSIONS AND SUMMARY 

- Antiadiabatic limit Our results can be compared with 
those obtained in the antiadiabatic limit. In fact, the 
antiadiabatic limit is relevant to rather few materials, but 
theoretically it provides a convenient benchmark. The 
antiadiabatic limit can be easily taken by integrating out 
phonons. The integration results in an effective retarded 
and attractive interaction between electrons. 



S, 



ee,eff 



= - f dTdT '^\ T y. 

J 2M 1 ^ 

x [n(r)-l][n(r')-l]. 



-iuj(t — t') 

j 2 + n 2 - 



(59) 



If the phonon frequency f2 is larger than typical elec- 
tronic energy scale such as A, then the frequency depen- 
dence in the denominator of e ^2+^2 - can be neglected. 
The summation over frequency gives the delta function 
in imaginary time, and Eq. (|59ll becomes 



-2E T 



drn T n x , 



n > a. 



(60) 



Now S ee ^ e ff can be combined with the on-site Coulomb 
repulsion U term, leading to 



U eS = U - 2E L 



(61) 



Note that this effective Anderson model is valid only in 
the limit f2 3> A. In our regime of interest U 3> El, U c s 
is sufficiently large so that the physics in Kondo regime 
applies. One can expect a qualitative change of behav- 
ior if U c s — U — 2El becomes negative. This is indeed 
confirmed by NRG studies 1^ Our treatment is not appli- 
cable to those regimes, and the polaron physics describes 
the regime adequately!^ The Kondo temperature of this 
effective Anderson model in the limit U S> El is given 



T K . 



off 



U cS A irU eS . 
— exph ^A- ] 



(62) 



The change of Kondo temperature due to electron- 
phonon interaction given by 



ST K . 



efi 



T, 



(0)/7T 

V4 



K 



A - 



A 



(63) 



in the limit of small A. Eq. (|63|l appears to be much 
larger than Eq. (|36II39|I . This also demonstrates that the 
physics in our regime is essentially different from that of 
antiadiabatic limit. 

-Adiabatic Expansion Deppeler and Millis studied the 
problem of electron-phonon interactions in correlated 
systems using adiabatic expansion of DMFT without any 
specific ground state assumed i&iii They have obtained 
the phonon Green function 



D- 1 (iu) = 1 



A 



lo 2 + \~fa p \uj\+o(~/ 2 ). (64) 
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The factor 1 — -A- represents the softening of phonon fre- 
quency in our notations. At half-filling and in the Hatree- 
Fock approximation, they obtained A c oc 1/u 3 while our 
result is A c oc 1/u . This is not real contradiction since 
our result does not apply to the half-filling case in the 
limit of large U . They have also identified the renormal- 
ized expansion parameters 



7 = 7 (1 -A/A,) 1 / 2 , A = 



A 



1 - A/A c 



(65) 



which coincide with our results. In our case, basically we 
obtain such expansion parameter through 



dui 1 - 

— {Y(-iu)Y(+iu)) N=0 = -Xj. (66) 

Z7T I 



associated with the diverse forms of couplings with elec- 
trons exhibit rather rich physics. 

-Summary We have studied AH model in the frame- 
work of semiclassical approximation, and interpreted the 
results in the light of DMFT. We obtained the analytic 
results for the Kondo temperature which are renormal- 
ized by weak electron-phonon interaction and for the 
phonon Green function which are valid at frequency 
higher than Kondo temperature. We have compared our 
results with the existing NRG results and other related 
works, and they are consistent with each other in their 
common regime of validity. 
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J drdr'Y^fl^G, 1 ^ -t)Mt) 



dr 



—4> 2 S - s (n T - ni) + gQ{ni - n L ) 



M /n N9 Am 2 , 
y(drQ) 2 + — 2~ <9 



(67) 



In this system the real spin is fully polarized, and the 
index a denotes the orbital degrees of freedom of ri- 
band |= d x 2_ y 2 1 {= e? 32 2„ r 2. Here Q is one of the 
Jahn- Teller phonons (the third component). In spite 
of the similarity of the local action Eq. (|67|l with the 
Eq. I|10(l except for the difference in the electron-phonon 
coupling, the outcomes are drastically different. For 
the model Eq. (|67|l . the phonon frequency is hardened 
wq = fly/ (U + AEl)/U, and Kondo temperature de- 
creases as follows: 



T K ~ Acxp[-i? 



U + 4E L 



], R > 1. 



(68) 



The reduction factor R is related to the overlap of the 
phonon wave functions in the process of instanton tun- 
neling. In fact, major results are just opposite to ours. 
This indicates that the orbital degrees of freedom being 
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APPENDIX A: THE DETAILED FORM OF Sxy 

In this section we exhibit all the terms comprising 
Sxy- The expansion of the action as a power series in Y 
is facilitated by utilizing the following parametric integral 



la 



o [1 + £C 2 (t)][1+£C 2 (t')] 



1 



In 



C 2 (r)-C 2 (r') 1 + C 2 (r') 



1 + Cg(r) 



(Al) 



Sxy consists of 10 terms which are numbered as 
SxY.i,i = 1, • • • , 10. Every term contains two time 
derivatives except for Sxy.i- 



A 



~>XY,1 



Y 2 (r) 



it J 1 + A 2 (r)' 



(A2) 



S 



XY,2 — 7 



drdr' \n\r -t'\Y(t)Y(t') 
(Jo ^ [ 1 + £A2(t)1 + £A2(t') ~ (l + W 1 



drdr' In |r - r'| Y(t)Y(t') ( / d£ 



2£X 2 (r) 



L (l+£A2(r)) 2 (l + £X2( T ')) (1+C% 2 ) 



2\3J 



(A3) 
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Sxy,3 = T / drdr' In \t - t'\X{t)X{t') 



1 1 

? [l + ^ 2 (r)] 3 1 ' [l + £A 2 (r') 

1 F 2 (t') (3A 2 (r')g 2 -g) 
? [l+£A 2 (r)] [l + £A 2 (r')] 3 

1 4^(r)X(^)r(r)y(rQ 
C [l+^ 2 (r)] 2 [l+^ 2 (r')] 2 

1 



(A4) 



?XY,4 = 



7T- 



drdr' ]h\t-t'\X(t')X(t) 



1 ^_d f(3X 2 (r)£ 2 -0A ^ 2 (r) 



+ / 

'o 



dr V [l + ^ 2 (r)]3 y [l+£X 2 (r')] 

1 Ap ± ( 1 \ (3X 2 (t')£ 2 -£)^V) 

"dr U + ^ a (r), 



[i + ^ 2 (t')] 3 

X(r) \ X(r')F(r)r(r') 



SxY,S = , 



dr V[l + £*2( r )] 2 / [l + £A 2 (r')] 2 
1 



A5) 



7T- 



drdr' ln|T-r'| A(t')F(t) 



X(r) 



F(r) 



X(r')Y(r') 



^dr V[l + ^ 2 (V.)] j ; [1 +a"-(V)] 

1 dearer 1 

°* ? dr ^1 + £X 2 (t)7 [1+£X 2 (t')] 2 



Sxvfi = j / drdr' \tl\t — t'\X(t') X(t) 



(A6) 



1 3X 2 (rK 2 -e^ 2 (r) 1 



[l + £A 2 (r)P dr l + £A 2 (r') 



d£4£ 2 



A(r)A(r') ^F(r') 
[1 + £A 2 (t)] 2 [1+^ 2 (t')F 



,(A7) 



Sxy,7 = j / drdr' In |r - r'| F(r') y(r) 



d 



1 



1 



5 



XY,S 



l + £A 2 (r)J l + ^ 2 (r')' 
' ' drdr' In |r - r' I 



1 



(A8) 



x[X(t>)Y(t)+Y(t')X(t)] 



[l+£A 2 (r)] 2 l + £A 2 (r') 
1 2£X(r')Y(r') 



+ X ^l + £Jf 2 (r) [l+£A 2 (r') 



5xf,9 = -j / drdr' In |r - r'| A(r')F(r) 



(A9) 



1 



[l+£A 2 (r)] 2 [l + £X 2 (r') 



(A10) 



Sxy,w = ^ I drdr' In |r - r'| F(r')X(r) 



*1 



1 



o [l+^ 2 (r')] 



2£A(r) 



o 



dr V[l + CA 2 (r)]; [l + £A 2 (r')] 



dr V[l+£A 2 (r)] 2 
2£X(r')Y(r') 



• (AH) 



APPENDIX B: DERIVATION OF RG 
EQUATIONS 



We will closely follow the original work by Anderson, 
Yuval, and Hamann^4, so that only the major steps will 
be described below. First the trajectories with instanton- 
antiinstanton pairs in the interval [to, to + dro] are to be 
integrated out. Let SZ be the part of partition function 
containing only one close pair. The trajectories contain- 
ing two or more pairs are accompanied with a factor of 
(dr ) 2 and are neglected in our (essentially one- loop) ap- 
proximation. 



SZ 



n[ 



l + z 



ti+l—To 
1 dt' 



t' -TQ 



dt" 



'*i+2r Jt>-[T +dTa] 

x exp[8V(t' - tx, ■ ■ ■ , t! - t")}] . (Bl) 



t',t" denote the locations of the close pair. Writing out 
the exponent V in detail 



SV 



K 1 J](-) l+1+J In 

j 

3 

K 2 J2 In 



t" - U 



TO 

t' - u 



TO 



t" - u 



K 2 J2 In 



TO 

a - 1 3 



TO 



K x In 



H - ?' 



d(t" - tj) 

t' - t" 



+ K 2 ln 



TO 



d(f - t"){B2) 



The last line of Eq. ()B2|) can be neglected since t — t" 
is the order of cutoff to. The most important difference 
between K\ part and K 2 part of Eq. (|B2|I is the absence 
of alternating sign factor in K 2 part. With the mean 
value theorem applied for the integration over t" of Eq. 
(|B1|) one can effectively impose the relation t" — t' — r 
in other parts of integrand, and then one factor of dro is 
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multiplied. Now the exponent V becomes 



sv = k^U 



i+.i 



In 



t' - 1, 



-In 



In 



In 



<' - ^ 



to 



ro 

d(t' - tj) 



f - 7& - tj 



ro 



t' —TQ — tj 



7ft 



dft'-TD-ij) 



(B3) 



We proceed by expanding e 5y ~ 1 + The Taylor 
expansion of the second term of K 2 part of Eq. I|B3|) 
respect to tq gives 

J2 J dt'd t ,d{t' - tj) = J2 [d(U+i - t s ) - d(ti - tj) 
3 3 
which vanishes upon the summation over i of the follow- 
ing approximation. 

JJ(1 + dr <h) ~ 1 + dT o Oi + o((dT ) 2 ). (B4) 

i i 

Collecting the remaining factors we obtain 

SZ = Yl [l + z 2 dT {(t t+1 -h-3T ) 

i 

rti+i—To 



'u+2t 



t' - ti 



2K 2 In 



t' - U 



TO 



d(t'-tj)]} . (B5) 



Re-exponentiating the expression in the big bracket of 
Eq. (|B5|I the correction to the partition function becomes 



SZ ~ exp ^z 2 dr (t i+ i - U - 3r ) 

i 



ti+1 — tj 



ti tj 



+2z 2 dr K 2 J2 dt ' ln 



t' - U 



d(t'-tff6) 



The first term in the exponent of Eq. IjBBjl contributes 
to the normalization of free energy e z /3dr ° . The K\ , K 2 
part of Eq. JB6|) behave very differently owing to the 
presence/absence of the oscillating factor. 

5Z ~ exp [z 2 dr ( J ft: 1 T )^(-) i+: ''(-l)21n|^ - tj 

i*3 



x exp 



dr { / dt'Y ( 2K 2 In 



t' - U 



tq 



d(t'-tj$7) 



The K 2 part of Eq. (|B7(I . in fact, does not depend on tj 
owing the integration over the whole range of imaginary 
time. Thus, the K 2 part of Eq. I|B7[) is proportional to 
the number of instantons, and consequently it contributes 
only to the fugacity correction. From the structure of Eq. 
(|B7ll the renormalization of K\ is easily found 



2 dr 



Ki = Ki — 4Ki(zt ) — - 
to 



(B8) 



The renormalization of fugacity is given by 



I ,1 \~ I ,T + dT 
(Tq + dT Q )z = (TqZ) X 



7ft 



exp 

where g 2 is 



Ki r + dr 
— — In 

2 t 



2 T dT (2K 2 g 2 ) , (B9) 



f^ln 


t' 






/ T 


TQ 



d(t'). (BIO) 

The dependence of d(t') is far stronger than ln ^- , thus 

we neglect the logarithmic factor up to logarithmic accu- 
racy. The explicit expression of d(t') is 



d(f) 



2A 7 2 A 



uj 2 e- luJt ' 



, j_ x 2tt u 2 + E*\u\ +u 2 ' 
Since t' ^ the above can be rewritten as 

2A 7 2 A r °° A ' ( TP*\i J 1 ,„2\„-iut' 



d(t') = 



^0 



duj (E*\uj\ + uj^)e- 
2^ uj 2 + E*\u\ +u)l 



(Bll) 



(B12) 



The is a small perturbation at both high and low 

energy. The correction coming from E* \u>\ is the order of 
E* /ujq <C 1. In leading approximation we can disregard 
E* . Explicit evaluation gives 



wo 



The evaluation of g 2 gives 

fj2 



^o T o 



(B13) 



(B14) 



The scaling equation of fugacity is 

^ = (zto) [(1 - ^) + (ztq) 2 (2K 292 )] . (B15) 
dr /to 2 

It is convenient to introduce new variables which are 
more suitable in the context of original Kondo model. 



dl — dro/ro, x = ztq, Ki = 2 — y, y <C 1. (B16) 
In terms of new variables the Eq. (|B15() can be written 



dx 
~dl 



y - 



x 2 ae~ UJoTo 



(B17) 



In the Eq. (|B17|) the short time cutoff was identified with 
1/A, while for the path integral approach the apparent 
cutoff is 1/U. This ambiguity has already been discussed 
in detail by Hamann in Sec. IV. B of Ref. ||, and it 
was argued that the ambiguity does not give rise to any 
significant errors. 
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